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ABSTRACT 

We critically examine the constraints imposed by carbon-enhanced metal- 
poor (CEMP) stars on the mixing mechanisms that operate in red giants. CEMP 
stars are created when the surface layers of a metal-poor dwarf are enriched with 
He-burning products via mass transfer from an evolved donor. The difference 
between main sequence (MS) and red giant CEMP abundances can be used as a 
diagnostic of the timescale for the mixing of the processed material into stellar 
interiors on the MS. Abundance trends with luminosity among red giant CEMP 
stars test theories of canonical extra mixing for low mass giants with a high 
bulk metallicity. We find a significant dilution in CN enrichment in giant CEMP 
stars relative to their MS precursors, and take this as evidence that thermoha- 
line mixing induced by mean molecular weight inversions is ineffective in CEMP 
stars. This contradicts models that rely on efficient thermohaline mixing in- 
duced by small fi gradients in red giants, because such models would predict that 
MS CEMP stars with large /i inversions would be homogenized on a very short 
timescale. The data do not rule out slower MS thermohaline mixing comparable 
to previously published estimates. We also find that canonical extra mixing is 
strongly suppressed in CEMP giants relative to stars with the same iron abun- 
dance. A likely cause is the shift in the location of non-equilibrium CN processing 
relative to the steep \i gradient in the hydrogen burning shell, which also occurs 
in solar metallicity RGB stars. Implications for the mass accreted by CEMP 
stars and the mechanism for canonical extra mixing are discussed. 
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Introduction 



It has long been known that asymptotic giant branch (AGB) stars can have their sur- 
faces enriched in helium, carbon a nd slow neutr on capture products, and such stars are 
known nucleosynthetic sources (e.g.. lHerwigll2005l ). If there is a nearby main sequence (MS) 



companion, mass transfer from the giant can produce surface abundance anomalies that 
persist even after the AGB star is gone. Such a system would appear to be a chemically pe- 
culiar MS star. In recent years there have been large-scale surveys for metal-poor stars that 
have included substantial spectroscopic follow-up campaigns, such as the Hamburg/ESO 
R-process Enhanced Star Survey (HERES: Christlieb et al. 2004). A large population of 
carbon-enhanced metal-po or (CEMP) stars has been discovered (e.g., iBeers fc Christlieb 



2005 



Lucatello et al. 120061) . and a majority of the CEMP stars have white dwarf compan- 



ions ( Lucatello et al. 20051 ) . This is strong evidence that there is a significant population of 



low mass metal-poor stars which have experienced mass transfer from a former AGB com- 



Frebel et al. 
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Tumlinson 
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(e.g., 

tant property of CEMP stars: they are unique astrophysical laboratories that can be used 
to test theories of mixing in low mass stars. A comparison of MS and evolved CEMP stars 
can test the timescale for mixing on the MS, and abundance trends within evolved stars 
can be used to test mixing on the red giant branch (RGB). In this paper we demonstrate 
that CEMP stars are not fully homogenized on the MS and that they do not experience the 
same RGB mixing pattern as other stars with the same iron abundance, and explore the 
theoretical explanations and consequences of both results. 

There are sound theoretical reasons to expect internal mixing in stars that accrete mass 
from an evolved companion. Metal-poor turnoff stars have thin surface convection zones, 
and as a result their surface abundances after mass transfer will reflect the composition of 
the accreted material. The He and C rich ejecta of an evolved AGB star has a higher mean 
molecular weight /i than a normal stellar envelope. Although the strong density gradients in 
a star prevent a classical Rayleigh- Taylor instability when such ejecta is accreted, Stancliffe 
et al. (2007) pointed out that this structure is unstable against thermohaline convection^. 
The precise timescale, however, depends critically on assumptions about the fluid geometry . 
Under the assumption of a simple spherical geometry (following iKippenhahn et al.lll980l ). 
Stancliffe et al. derived a characteristic timescale of order 1 Gyr for homogenizing a CEMP 
star. In this theoretical model the surface abundances are diluted through 90% of the mass 



1 Thermohaline convection is driven by a density difference between a fluid clement and its surround- 
ing medium, provided that the latter has a higher mean molecular weight and the heat exchange keeps a 
temperature difference sufficiently small. 
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of the star within this timescale for typical companion masses and ages, with only the 
nuclear processed core escaping full mixing. As supporting evidence the authors cite both 
the absence of extremely high ([C/Fej^|> 2.5) enrichment in MS stars and the apparent lack 
of a difference between MS and RGB surface abundances. 

There has also been considerable interest in thermohaline mixing as an explanation 
for deep mixing in the envelopes of red giants. There is clear evidence that ordinary 
low-mass giants experience in situ mixing of CN- processed material as they become more 
luminous (e.g.. iGratton et al.l 120001 ; ISmith fc Martelll 120031 ) . Because of its universality, 



Denissenkov fc VandenBer 3 63) proposed to call it canonical extra mixing. Rotational 



mixing has long been a popular explanation fe.g.JSweigart & Mengellll979uCharbonne 



19981 ; iDenissenkov &: Toutll2000l ; iDenissenkov &: VandenB erg 
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Denissenkov et al 



Brown, fc Wallerst 
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but v i able models have pro ven to be difficult to construct (jChaname. Pinsonneault. &: Terndrup 
20051 ; iPalacios et al.ll2006l ). Attention has therefore recently shifted towards other physical 
mechanisms. 

When performing hydrodynamic sim ulations using the 3D stellar evolution code D j ehuty, 
Eggleton. Dearborn, fc Lattanzid (120061 ) found unexpected convective-like motions in the ra- 
diative envelope of an upper RGB model. They attributed this to a Rayleigh- Taylor insta- 
bility (RTI) driven by a mild (A/i//i « — 10 -4 ) /J inversion generated by 3 He burning. 3 He 
is produced by non-equilibrium pp burning on the MS and dredged up on the lower RGB. 
On the upper RGB material enriched in 3 He falls onto the hydrogen burning shell and the 
3 He( 3 He, 2p) 4 He reaction lowers /x, while equilibrium H burning raises \i in still deeper layers; 
the net effect is a well-defined local minimum above the hydrogen burning shell. 

We note that the RTI would not be expected in compressible, and hence stratified, 
stellar material. Instead, thermohaline convection (for the same reasons as in the CEMP 
case) would be a natural consequence. This phenomenon would be an attractive candi- 
date to drive canonical extra mixing if it could provi de mixing depths and rates consistent 
with their empirically constr a ined values. However, ICharbonnel &: Zahnl (120071 ) as well as 
Denissenkov &: Pinsonneaultl (120071 ) have demonstrated that the 3 He-driven thermohaline 
convection will be fast enough to explain the RGB mixing patterns only if fluid elements 
have finger- like structures with a length to diameter ratio l/d > 10. This produces a much 
shorter mixing timescale than proposed by Kippenhahn et al. (1980). 

If canonical extra mixing is to be explained by thermohaline convection, the inferred 
mixing timescale in MS CEMP stars should be extremely short (of order 1-10 Myr). If 



2 We use the standard spectroscopic notation: [A/B] =log (JV(A)/JV(B))-log (iV(A)/JV(B)) Q , where N(A) 
and N(B) are number densities of the nuclides A and B. 
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CEMP stars are fully mixed after a mass transfer episode, this implies that there will be no 
dilution of the abundance anomalies during the classical first dredge-up (FDU). By contrast, 
if mixing is inefficient on the MS one might expect that large abundance anomalies confined 
to the outer layers will be diluted when the star goes to the RGB and develops a deep 
surface convection zone. We therefore propose to use the relative abundances in MS and 
RGB CEMP stars as a diagnostic of the mixing timescale on the MS, and by extension as a 
test of the timescale for thermohaline convection in a very favorable environment. 

Canonical extra mixing might also be expected to set in for CEMP stars when they reach 
the upper RGB, and a distinctive manifestation of canonical extra mixing at low metallicities 
is a decline of [C/ Fe] with increased luminosity (e.g. Sm ith & Martel 2003). There have been 
published claims ( ILucatello et al.ll2006l ; lAoki at al.ll2007l ) that the abundances of CEMP stars 
decline with increased luminosity. However, these studies did not separate out the MS/RGB 
transition (which tests the depth of accreted material) from the RGB trends. We therefore 
both examine the data for trends with luminosity and compare to theoretical expectations 
for empirical models based on normal metal-poor giant data. Again, in addition to the direct 
implications our test can be used to study the metallicity dependence of RGB mixing. It is 
known that solar-ab undance open-cluster RGB stars have weake r extra mixing than globular- 



cluster giants (e.g., Chaname. Pinsonneault. fc Terndrupl bood ) . However, the open-cluster 



giants are both more metal-rich and more massive; both properties could influence the mixing 
pattern. CEMP stars permit a direct comparison of low-mass stars with different metal 
abundances. 



2. Observational Data 

Ideally, our goal would be to select a sample of binary mass transfer systems where 
we can look for abundance trends both between MS and RGB stars and with luminosity 
on the RGB. There are now large samples of CEMP stars, but there are obvious poten- 
tial biases associated with examining abundance patterns in samples defined by abundance 
measurements. As a result, some discussion of the selection of the observational sample is 
essential. 

One important clue is that the amount of carbon in CEMP stars is very high relative 
to halo non-C-enhanced sta rs. This is most heavily seen in the anti-correlation between 



[C/Fe] and [Fe/H] (Fig. 3 in lAoki at al.ll2007l ). This is the pattern that one would expect if 



mass transfer from an evolved companion gave the recipient star so much C that the initial 
inventory was negligible. However, especially at the iron-poor end, some chemical evolution 
models might be able to produce exotic abundance patterns for the primordial composition 
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(e.g.. lRyan et al.ll2005l ; iFrebel et al.ll2005l : iKomiya et al.ll2007l ). Fortunately, the binary mass 
transfer model makes other predictions: there should still be a white dwarf companion and 
there should be associated s-process anomalies. A majority of CEMP stars exhibit radial 
velocity variations consistent with a white dwarf companion; the observed fraction of 68% is 
consistent with all CEMP stars being in bi naries (Lucatello et al. 2005). Extensiv e surveys 
have also shown s-process enrichment (e.g., iBarklem at al.ll2005l ; lAoki at al.l 120071 ). 



We have therefore used two samples of CEMP stars: a smaller sample from lAoki at al. 



( 120071 ) in wh ich only stars with s-p rocess enhancement ([Ba/Fe] > +0.5) were included, and 
a larger one (jLucateUo et al.ll2006l ) where the standard definition ([C/Fe] > 1) to distinguish 
between CEMP and normal field halo stars was used. Lucatello et al.'s sample is in fact a 
subsample of stars from the HERES survey (IBarklem at al.ll2005l ) that were re-analyzed. We 
did not use the stars with [C/Fe] > 1 from Barklem et al.'s original paper. The automated 
procedure used in that paper to infer T e g values from colors became unreliable for CEMP 
stars with strong CH, C2, and CN molecular bands (Lucatello, 2007, private communication). 
The luminosities were inferred from the surface gravity and T e g with an assumed mass of 
M = 0.8-0.85 M Q . The values of log L/L & < 0.8 and log L/L & > 1.4 were defined as the 
divisions between unevolved MS and subgiant stars that have not experienced the FDU yet 
and evolved RGB stars in which the FDU has already ended. We have found that 10 out 
of 14 unevolved stars in the second sample have [Ba/Fe] > +1 according to Barklem et al.'s 
data, which is a distinctive feature of accretion of AGB material. 

Under the binary mass transfer hypothesis the initial carbon inventory is dwarfed by 
the transferred carbon, so the physical quantities of interest are [C/H] and [N/H] rather than 
[C/Fe] and [N/Fe]. As a result, we work in the former plane. There is a striking difference 
in carbon as a function of luminosity between CEMP and normal halo stars (Fig. [TJ blue 
open and black filled circles, respectively), with canonical extra mixing being clearly seen 
in the non-C-enhanced stars above log L/L® pa 2.2 but not in the CEMP stars. The black 
curve reproduces the evolutionary decline of [C/H] in the non-C-enhanced stars by our stellar 
models with canonical extra mixing (for details, see next section). The red curve is a model 
constructed under the assumption that neither the MS thermohaline convection nor RGB 
extra mixing have taken place. The only mixing event experienced by this model is the 
standard FDU. It seems to adequately describe the [C/H] abundances in the CEMP stars. 

We compare the trends in our small and large CEMP samples in Fig. [2J and find that 
they are very similar. In fact, the [C/H] abund ances in the two samples have the same 
linear fit (jLucateUo et al.ll2006l ; lAoki at al.ll2007l ). We also note that there is a small, but 



real, d ifference between the MS and RGB abundances. To be more precise, iLucatello et al. 



( 120061 ) have obtained very similar slopes -0.18 ± 0.05, -0.21 ± 0.05, and -0.19 ± 0.04 for 
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[C/H], [N/H], and [(C+N)/H], respectively. We argue that this is evidence for dilution of 
enriched material in the MS stars by the classical FDU in the RGB stars. Both effects are 
quantified below. 



3. Reduced Efficiency of Canonical Extra Mixing in RGB CEMP Stars 

fn this section, we show that the efficiency of canonical extra mixing in CEMP stars 
is strongly reduced. This is most likely due to shallower mixing which does not penetrate 
deep enough to dredge up carbon-deficient material. The most plausible interpretation of 
the shallower mixing depth i s that it is limited by a st eeper /i-gradient built up by H burning 
on the abundant 12 C nuclei. iDenissenkov et al.l (120061 ) came to a similar conclusion for solar- 
metallicity red giants. Given that RGB mixing is inefficient in CEMP stars, we can use the 
C and N abundance distributions for unevolved (the MS and subgiant branch) and evolved 
(upper RGB) CEMP stars to test the FDU and draw conclusions about the efficiency of 
thermohaline convection in MS CEMP stars. This will be done in the next section. 

The slope of [C/H] as a function of log L/L & produced by canonical extra mixing de- 
pends on the mixing depth and diffusion coefficient D m ix- We model canonical extra mixing 
by setting D mix equal to a fixed fraction of the thermal diffusivity K. This parameterization 
takes into account that large-scale mixing instabilities in a stably stratified radiative zone 
may only develop on a timescale longer than the thermal one, in which case one would 
expect -D m i x < K. The known examples of such instabilities that have been proposed 
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Charbonnel & Zahnl 120071; IDenissenkov & Pinsonneaultl 120071) 



the latter case, the diffusion coefficient can be estimated using the following heuristic ex- 
pression: 



D ■ ps 

^mix ~ 



IV, 



V a ,H - V 



■K, 



rati 



where V^, V a d, and V ra d are the /x- and temperature (adiabatic and radiative) gradients 
(logarithmic and with respect to pressure). Equation (ED) adequately represents b oth the 
diffusion coefficient for thermohaline convection derived by iKippenhahn et al.l (119801 ) 
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who argued that the length (or, the mean free path) I of a convective element should be of 
the order of its diameter d, and the rate of mixing by elongated narrow "fingers" (I > d) 
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advocated by lUlrichl (Il972l ) 



8 2 K 

-Duirich — -7T — — 

<J Vad — Vrad 
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(3) 



Unfortunately, the ratio l/d cannot be estimated from first princi ples. Semi-empirical 
models require lid > 10 to generate suffic ient mixing to explain the data (jCharbonnel &: Zahn 
20071 : iDenissenkov fc Pinsonneaultl 120071 ). The maximum depth of the 3 He-driven thermoha- 
line convection is uncertain also. The naive assumption that it co incides with the location of 
the m inimum on the yU-profile produces mixing that is too shallow (IDenissenkov &: Pinsonneault 
20071 ). Therefore, we have chosen to specify the mixing depth r mix in units of R Q and adjust 
it appropriately. This choice is motivated by our finding that the fi(r) does not change much 
with time. In the absence of any information on how the mixing depth should vary with 
time, keeping r mix constant is a reasonable parameterization. 

To allow an interpolation in the abundances of He and C as the AGB material gets 
redistributed inside an accreting star by convec tion and extra mixing we have modified our 



stellar evolution code (Denissenkov et al. 20061) to include the 2005 version of the OPAL 



equation of state ( Rogers. Swenson. fc Iglesiasl 1996)^1. Taking into ac count that the average 



metallicity of the non-C-enhanced stars from the spectroscopic study of lBarklem at al.l (120051 ) 
is ([Fe/H]) « —2.7, we have represented them with a stellar model having the heavy-element 
mass fraction Z = 4 x 10~ 5 . For the initial helium mass fraction, stellar mass, and mixing 
length parameter we have used the values Y = 0.24, M = 0.85 M , and «mlt = 1-58, 
respectively. 

We begin with the normal halo star sample, and we use them to constrain our model of 
canonical extra mixing. The sharp decline of [C/H] at log L/L Q > 2.2 in the non-C-enhanced 
very metal-poor (VMP) stars (Fig. [I]) has been reproduced by us with r mix = 0.05 R Q 
and -Dmix = 0.1 K. The chosen depth is shown with the vertical black line in Fig. HA. It 
guarantees that material deficient in C is dredged up to the bottom of convective envelope, 
as indicated by observations. A similar result can be obtained for a slightly smaller r mix with 
an appropriately reduced value of D mix , e.g. for r mix = 0.045 R® and D m[x = 0.04 K. At the 
adjusted mixing depths T(r mix ) is lower than the temperature at the base of the H shell by 
AlogT 0.24 and 0.19, respective ly, which agrees very well with the estimates obtained by 
Denissenkov fc VandenBergi (120031 ) based on globular-cluster data. 



Our results for the CEMP stars are plotted in Fig. [31 The short-dashed red and solid blue 
curves have been computed for a stellar model with Y = 0.24 and Z = 10~ 4 (the CEMP 



3 The update is at http://www-phys.llnl.gov/Research/0PAL/EQS_2005/ 
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stars have ([Fe/H]) ~ —2.1) whose initial mass and CN abundan ces were M = 0.7M e 



[C/Fe] = +0.3, and [N/Fe] = -0.5. Following IStancliffe et all (120071 ). we have assumed that 



0) 



at the age of 0.86 Gyr the star began accreting material with the rate 10 6 M yr 1 from 
its AGB binary companion. The total mass accreted and its He and CN abundances were 
chosen to be M accr = 0.15 M , Y = 0.34, [C/Fe] accr = +2.1, and [N/Fe] accr = +1.6. 

Canonical extra mixing in our CEMP model has been modeled using the same mixing 
parameters that we have adjusted for the non-C-enhanced VMP stars. Comparing the [C/H] 
data at log L/L & > 2.4 with the short-dashed red and solid blue curves in Fig. [3]A., we infer 
that the efficiency of canonical extra mixing appears to be strongly reduced in the CEMP 
stars compared to their non-C-enhanced VMP counterparts. 

The reduced efficiency of canonical extra mixing in the CEMP stars becomes even more 
evident when we consider their [N/H] values as a function of luminosity (Fig. EfB). CN 
processing for a high initial carbon should result in a sharp increase of [N/H] on the upper 
RGB (the short-dashed red and solid blue curves at log L/L Q > 2.4). Contrary to this, the 
observed trend of [N/ H], as well as th o se of [C/H] and the sum [(C+N)/H], with luminosity 
were characterized by Lucatello et al. ( 20061 ) as gradual declines (long-dashed black lines in 



Fig. [3]). They have ascribed them to a dilution "by mixing with a much larger H content". 
We concur in part, and interpret the decline from the MS to the RGB as being a dilution 
effect produced by the FDU alone. However, we also believe that representing this dilution 
as a continuing process occurring on the RGB is not supported. The FDU is effectively 
complete at log L/L Q ~ 1.4. At 1.4 < log L/L Q < 1.8 the convective envelope mass is still 
increasing, but only by another 3-4%. So, no noticeable changes of the surface chemical 
composition due to the FDU dilution would be expected above log L/L Q « 1.4, and we 
contend that no trend is seen. 

Canonical extra mixing begins to operate only at log L/L Q rs 2.4 but it must increase 
the surface N abundance, as it really does in non-C-enhanced VMP stars (the open blue 
circles in Fig. [5]), and not decrease it. Hence, a dilution of the accreted He-, C-, and N-rich 
material can only occur either on the MS, as a result of thermohaline convection, or on the 
lower RGB at 0.8 < log L/L Q < 1.4, as a result of the FDU. We will return to this problem 
in the next section. Here, we have to conclude that, as revealed by the non-increase of [N/H] 
in the upper RGB CEMP stars, canonical extra mixing does not manifest itself in them. This 
conclusion agrees with the different behavior of [C/H] at high luminosities in the CEMP and 
non-C-enhanced VMP stars (Fig. [1]). 

Interestingly, there is another group of stars in which canonical extra mixing is observed 
to have a low efficiency. These are low-mass RGB stars of nearly solar metallicity. Only mod- 
est depletions of the carbon isotopic ratio ( 12 C/ 13 C ~ 10-15) are reported in them, while 
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the variations of [C /H] and [N/H] are close to their FDU limits (e.g.. iGilroy &: Brownlll991 



Origlia et al.ll2006l ). This contrasts with metal-poor red giants in which canonical extra mix- 



ing often reduces the 12 C/ 13 C ratio down to its e quilibrium value ~4, which is accompanied 



by the well-marked evolutionary decline of [C/H]. ICharbonnel. Brown, fc Wallersteinl (119981 ) 
have shown that such dependence of the efficiency of canonical extra mixing on metallicity 
is expecte d if the mixing cannot pe netrate below some universal critical //-gradient (V At ) cr it- 



Recently, iDenissenkov et aJj (120061 ) have proposed that the shallower mixing in the solar- 
metallicity red giants is caused by a contribution to the //-gradient from the H depletion 
that accompanies the C to N transformation. This contribution increases with the metal- 
licity because N(C) oc Z. At a metallicity close to the solar one, the H burning on the 
12 C nuclei builds up an additional hump in the /i-gradient profile that prevents mixing from 
approaching the H-burning shell as close as in metal-poor giants. If this interpretation is 
correct then the reduced efficiency of canonical extra mixing in RGB CEMP stars is almost 
definitely caused by their much higher C abundances compared to their non-C-enhanced 
counterparts. 

To test this hypothesis, we confront the CN-abundance and profiles in the vicinity 
of the H-burning shell in our two upper RGB models in which the H shell has just erased the 
discontinuity in the H-abundance profile left behind by the bottom of convective envelope 
at the end of the FDU (it is not until this moment when canonical extra mixing begins to 
manifest itself). The first model represents the non-C-enhanced VMP stars: it has M = 
0.85 M Q , Z = 4x 10~ 5 , [C/Fe] = +0.3, and [N/Fe] = -0.5 (panels A and C in Fig. HJ. The 
second had the initial mass M = 0.7 M and abundances Z = 10~ 4 , [C/Fe] = +0.3, and 
[N/Fe] = -0.5 but it had accreted 0.15 M© of material with [C/Fe] = +2.1 and [N/Fe] = +1.6 
(panels B and D in Fig. H]). The vertical line crossing panels A and C shows the mixing 
depth r mix = 0.05 that has been used to model the evolutionary decline of [C/H] in the 
non-C-enhanced VMP stars. The horizontal line crossing panels C and D reads the value of 
(V M ) C rit ~ 3.6 x 10 -5 corresponding to this mixing depth. When applied to the C-enhanced 
model this (V M ) CT it determines a shallower mixing depth (r m i x « 0.061, panels B and D) with 
which extra mixing will only be able to diminish the surface 12 C/ 13 C ratio without changing 
[C/H] and [N/H] (panel B). For comparison, dotted curves in panels B and D show the C- 
abundance and V M profiles in a solar-metallicity 1 M RGB model. It is seen that in terms 
of the depth of canonical extra mixing the CEMP model is equivalent to the solar-metallicity 
model. Hence, it is most likely the mixing depth that is affected by the variation of the star's 
C content in the first place. 

For our CEMP model we have intentionally chosen the large positive value of [N/Fejaccr 
because it resembles the observed N abundance anomaly in unevolved CEMP stars. It is 
thought that 12 C dredged- up from the He-burning shell in AGB stars might have partially 
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been converted into N in the hot-bott om burning (in more massive AGB s t ars) or by extra 



2006 



Sivarani et al. 2006; Johnson et al 



mixing similar to the canonical one (Nollett. Busso. fc Wasserbura 120031 ; iMasseron et al. 



200 



j). To illu strate this poin t , we h ave reproduced 



in a slightly modified form Fig. 4 from the paper of ISivarani et al. (hood) in our F ig. 
For non-C-enhanced VMP stars, we have used the data obtained by lSpite et al.l (120051 ) (the 
filled and open blue circles denote unmixed and mixed RGB star s in th eir terminology) . For 
CEMP stars, we have used the data compiled by ISivarani et al.l (120061 ) (the filled and open 
red squares represent stars with log g > 3.5 and log g < 3.5, respectively). We have also 
plotted the C and N a v erage abundances returned to the ISM by metal-poor AGB stars, 
as predicted by iHerwigl (120041 ) (the magenta asterisks with numbers below them giving the 
AGB model's mass). The black curve in Fig. [5] shows changes of the surface C and N 
abundances in an RGB model with M = 0.83 M and Z = 0.0001 produced by canonical 
extra mixing with the depth r mix = 0.045 R Q and rate D mix = 0.04 K . It is important t o note 
that open blue circles {mixed stars) are located on the upper RGB (ISpite et al.ll2005l ). This 
observational fact and a comparison of ou r black curve with the blue data points support 
the conclusion made by ISpite et al.l (120051 ) that they had revealed the abundance patterns 
produced by canonical extra mixing in VMP stars. 

It is interesting that most of the CEMP stars in Fig. [5] fill the gap between the C and 
N abundance yields predicted for the low- and intermediate-mass AGB stars. The latter are 
well separated from the former because in the intermediate-mass AGB models C is efficiently 
transformed into N at the bottoms of their convective envelopes (the hot-bottom burning). 
The mysterious location of the CEMP stars in the gap probably indicates that in their low- 
mass AGB donors C was also transported down to a high-temperature region where it was 
partially transformed into N. The rathe r low carbon isot opic ratios ( 12 C/ 13 C< 40) in low- 
luminosity s-element-rich CEMP stars ( jRyan et al.l 120051 ) support this hypothesis because 
witho ut such extra mixing the low-mass AGB stars would produce 12 C/ 13 C > 1000 (IHerwig 
20041 ). 



4. Low Efficiency of Thermohaline Convection in MS CEMP Stars 

In the previous section, we have presented observational and theoretical arguments for 
the shallower depth of canonical extra mixing in RGB CEMP stars. The depth seems to 
be shallow enough for extra mixing to leave the post-FDU C and N abundances unchanged. 
As we noted, the FDU starts at log L/L Q « 0.8 and it effectively ends at log L/L e « 1.4 
(the short-dashed red and solid blue curves in Fig. [3]). Hence, we can compare the C and N 
abundance distributions for the CEMP stars with log Lj L® < 0.8 and for those of them with 
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log L/L & > 1.4 to estimate the degree of dilution of the CN-rich material by the FDU. These 
are plotted in panels A, B (for our larger sample), and C (for the smaller sample) in Fig. [6j 
We see that the distributions for the stars that have not experienced the FDU yet (shaded 
histograms) are shifted toward higher abundances with respect to the distributions for the 
stars in which the FDU has already ended (thick solid histograms). This holds true for both 
the C (panels A, C) and N (panel B) abundance. There is a difference in the mean too, but 
it is not statistically significant because of the large star-to-star scatter. However, there is 
a clearly defined median, and we employ a median based test. We interpret the abundance 
distribution shifts as a signature of the FDU dilution. The amplitude of the shift depends 
on both the average accreted mass (M accr ) and the maximum mass of convective envelope 
attained during the FDU M™^ x « 0.5 M . To estimate their ratio, and hence (M accr ), we 
shifted the shaded histograms to the left by small steps in [X/H], where X stands for C and 
N (the shifts were applied to each star in the group), and we calculated the Kolmogorov- 
Smirnov probability Pks that the shifted shaded and thick solid histograms did not differ 
from one another after each step. Results of this statistical test are plotted in Fig. 03. They 
show that the same shift A[X/H] ~ —0.4 leads to the maximum and quite large values of 
Pks f° r both C and N. Neglecting the initial CN abundances, as compared to their accreted 
and post-FDU values, we estimate (M accr ) « 10 A[X/H] Mgg* w 0.2 M Q . 

The fact that the pre-FDU and post-FDU C and N abundance distributions are surpris- 
ingly well separated implies a low efficiency of thermohaline convection in MS CEMP stars, 
which has immediate implications for RGB mixing. In Fig. [71 we illustrate the /^-profiles in 
our MS (red curve) and upper RGB (blue curve) CEMP models. On the scale of the MS 
inversion (A/x//x)ms ~ 0.1 the tiny RGB depression (A/x//x)rgb ~ 10~ 4 at r m 0.067 P Q can 
hardly be seen. When comparing these profiles, one should keep in mind that the efficiency 
of thermohaline convection is also proportional to the thermal diffusivity K (equation [1]). In 
the MS model, K is three to four orders of magnitude lower than in the radiative zone of the 
RGB model. As a result, the MS thermohaline diffusion coefficient can be less than 10% of 
-D m ix in the RGB model. However, this rate is still fast enough to mix the MS CEMP star 
on a timescale shorter than 10 Myr, i.e. almost instantaneously compared to the star's MS 
lifetime. The evidence that this does not happen casts doubt on thermohaline convection as 
the mechanism for canonical extra mixing. 

5. Conclusion 

In this paper, we have put forward the hypothesis that the anomalously high C abun- 
dances in the material accreted by MS CEMP stars result in shallower depths of canonical 
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extra mixing in these stars during their subsequent RGB evolution. The shallower mixing 
depths are limited by steeper /^-gradients produced by the H burning on abundant 12 C nuclei. 
A similar phenomenon occurs in open cluster RGB stars, which are of comparable total metal 
abundance. Until now, it was not clear whether the reduced efficiency of mixing in open 
cluster giants was caused by their higher metallicity or their higher mass relative to globular 
cluster giants. We now see the same trend in low mass CEMP giants, and can therefore con- 
clude that metallicity is the primary determinant of the depth of mixing. Mixing would still 
be expected in the outer envelopes of these giants, and could impact other diagnostics, such 
as the 12 C/ 13 C ratio. Observational data on abundance trends with luminosity of tracers of 
shallower mixing, such as the carbon isotopic ratio, would be a useful further test of this 
theoretical prediction. 

For thermohaline mixing to be effective on th e RGB the timescale must be 100 to 1000 



times shorter than the iKippenhahn et al. (1980) e stimat es. A comparable increase in ef- 



ficiency on the MS relative to the IStancliffe et al.l (120071 ) estimates would imply a mixing 



timescale far shorter than the MS lifetime. As a result, we would expect to see no dilution 
effect in the abundances, which contradicts the data. We therefore conclude that thermo- 
haline mixing is inefficient in MS CEMP stars, contrary to the natural prediction of models 
where it is the dominant mechanism for RGB mixing. In light of the potential importance of 
this conclusion, we now discuss the limits that we can place on the timescale of MS mixing, 
potential observational and theoretical issues, and future observational tests. 

We cannot make a simil arly strong claim for th ermohaline convection operating over Gyr 
timescales, as advocated by IStancliffe et al.l (120071 ). Our globally averaged result is a 0.2 M Q 



average depth of the accreted material; this could be an unmixed 0.2 M accreted layer, or 
an initially lower mass accreted layer that has diffused to an average depth of 0.2 M & . A 
model with a diffusive timescale of order the post-mass transfer MS lifetime and a 0.1 M & 
accreted layer, for example, would be consistent with our results. When the mixing timescale 
is comparable to the evolutionary lifetime, a population synthesis approach is required for 
a comprehensive analysis. One must account for the distribution of abundance anomalies 
in the donor stars, mass accretion amounts, and the range of secondary masses over which 
the MS recipient can become a red giant. Such calculations would be interesting, but are 
beyond the scope of this paper. 

Observational selection effects are also a concern. If there was a differential bias in the 
effective temperature estimates for dwarfs and giants, for example, the relative abundance 
patterns could be impacted. Similarly, a selection criterion where dwarfs would require 
higher carbon anomalies than giants to be discovered would manifest itself as a dilution 
effect. There is no evidence for the latter effect in the data; the observed carbon anomalies 
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are well-separated from the normal halo star pattern. Effective temperature scale estimates 
are a more significant concern, especially for color-based es timates, because o f the spectral 
energy distribution distortions from very high [C/Fe] ratios. lAoki at al.l (120071 ). for example, 
found evidence for iron abundances dependent upon the excitation potential of the lines, 
which is usually a signature of an incorrect effective temperature. However, to impact our 
results there would need to be a differential impact on the abundances derived in dwarfs and 
giants rather than an overall systematic effective temperature shift, and Aoki et al. did not 
report such an effect. We also note that very similar relative trends appear in samples where 
spectroscopic and color-based effective temperature estimates are employed. 

Is there a theoretical mechanis m that could differentially suppress thermohaline convec- 



tion in MS (but not RGB) stars? iDenissenkov &: Pinsonneault 



(120071 ) have proposed that 



thermohaline convection can be suppressed by horizontal turbulent diffusion. It is, at least 
in principle, possible to invoke stronger horizontal turbulence in dwarfs than in giants. An- 
gular momentum constraints from the existence of rapid rotation in horizontal branch stars 
do imply strong differential rotation with depth in the connective e nvelopes of RGB stars 
( jPetersonl Il983l : iPinsonneault et al.l Il992l ; ISills &: Pinsonneaultl 120001 ) . One would therefore 
expect rotation rates in the RGB envelope comparable to those in MS radiative interiors, 
and as a result comparable degrees of horizontal turbulence. This is an area that deserves 
further theoretical scrutiny. 

There is indirect evidence for extra mixing in low-mass AGB stars (e.g., see Fig. 
and its discussion in section 3). However, if extra mixing in RGB stars is really driven by 
3 He burning then it should die out by the end of the RGB evolution because of the 3 He 
exhaustion. In this case, the 3 He-driven thermohaline convection could not resume working 
in the same stars on the asymptotic giant branch (AGB). So, we would expect the absence 
of observational signatures of extra mixing in low-mass (M < 2M ) AGB stars unless the 
mixing in them is of a different nature. However, given the similarities in their depth and 
in the structure of radiative zone where they operate, it is unlikely that the RGB and AGB 
mixing have different physical mechanisms. 

Further observational work would be extremely fruitful. Samples should be examined 
for their binary nature in order to generate binary mass transfer samples, and both the 
reliability of temperature estimates and the relative and absolute abundance errors should 
be rigorously examined. If thermohaline convection occurs over an intermediate timescale, 
there should be a difference in the mean abundance level between MS stars that experienced 
mass transfer recently and stars with older mass transfer episodes. Hotter white dwarfs can 
be detected photometrically, particularly for dwarf CEMP stars, and a difference in sur- 
face abundance between systems with younger and older white dwarf companions would be 
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definitive evidence for (or argument against) a moderate timescale for thermohaline mixing. 

We are grateful to Sara Lucatello for sending us the C and N abundance data for 
the CEMP stars. It is our pleasure to thank Jennifer Johnson, Thomas Masseron, Grant 
Newsham, Kristen Sellgren, and Don Terndrup for usefull discussions. We acknowledge 
support from the NASA Grant No. NNG05 GG20G. 
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Fig. 1. — The C abundance as a function of luminosity in the CEMP (black filled circles) 
and non-C-enhanced VMP (blue open circles; the blue cross shows the errorbars) stars. The 
latter are seen to experience canonical extra mixing at log L/L Q > 2.2. The declining part of 
the black curve has been computed for the mixing depth r mix = 0.05 R Q (panels A and C in 
Fig. H]) and diffusion coefficient D mix = 0.1 K, where K is the thermal diffusivity. The only 
mixing episode experienced by the CEMP stars appears to be the standard first dredge-up 
at 0.8 < logL/L© < 1.4 (red curve). 
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Fig. 2 . — The [C/H] as a function of luminosity for Ba-enhanced CEM P stars fromlAo 



(120071 ) (open circles) and for a larger sample of CEMP stars from iLucatello et al 



d at al. 
J2006h 



(filled circles). The cross shows the data errorbars. Both samples have the same linear fit. 
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Fig. 3. — The C and N abundances in the CEMP stars (circles) are compared with those 
predicted by our stellar evolution models in which canonical extra mixing has been pa- 
rameterized with r mix = 0.05, D mix = 0.1K (short-dashed red curves), and r mix = 0.045, 
-Dmix = 0.04 K (solid blue curves). Crosses show the observational errorbars. Black long- 
dashed lines are linear fits to the data with slope s — 0.18±0.05, — 0.21± 0.05, and — 0.19±0.04 
for [C/H], [N/H], and [(C+N)/H], respectively (Licatello et alJboOfih . 
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Fig. 4. — The abundance (panels A and B) and /^-gradient (panels C and D) profiles in our 
non-C-enhanced (panels A and C) and C-enhanced (panels B and D) upper RGB models. 
Dotted curves in panels B and D are from our solar-metallicity 1 M® upper RGB model. 
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Fig. 5. — The C and N abu ndances for the u nmixed and mixed VMP stars (filled and 
open blue circles) obtained by ISpite et al.l (120051 ) as well as those for the CEMP stars wit h 
log g > 3.5 and log g < 3.5 (filled and open red squares) compiled by ISivarani et al.l (120061) . 
Magenta asterisks represent the average abundance yields predicted by iHerwid (12004 ) for 
metal-poor AGB models whose masses (in M©) are given below the data points. Black curve 
shows changes of the surface C and N abundances produced by canonical extra mixing with 
the depth r mix = 0.045 R Q and rate D mix = 0.04 K in an RGB model with M = 0.83 M 
and Z = 0.0001. Note that most of the CEMP stars are located in the gap between the 
predicted AGB yields. 
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Fig. 6. — The C (panels A and C) and N (panel B) abundance distributions for the 
CEMP stars that have not experienced the FDU yet (shaded histograms) and for those 
of them in wh i ch the FDU has already ended (th ick solid curves). The data are from 
Lucatello et al.l (120061 ) (panels A and B) and from lAoki at al.l (120071 ) (panel C). Panel D 
shows the Kolmogorov-Smirnov probabilities that the shaded histograms in panels A and B 
do not differ from their thick solid counterparts after being shifted to the left by A [X/H]. 
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Fig. 7. — The /x profiles in two models with M = 0.85 M Q and Z = 0.0001. The first 
(red curve) is our CEMP MS model that began its evolution as a 0.7 M star and has just 
accreted 0.15 M & of AGB material enriched in He and C. The second (blue curve) is a model 
in which the H shell has just erased the H-abundance discontinuity. It has preserved its 
initial composition. The blue curve has a local minimum of Afi/fi « 1.4 x 10" 4 located 
at r ~ 0.067 -Rq. It is produced by 3 He burning. If canonical extra mixing in upper RGB 
stars is identified with thermohaline convection driven by this tiny /i inversion then a similar 
mechanism should produce very efficient mixing in CEMP MS stars as well, given their 
much stronger /i inversions (red curve), even though the thermal diffusivity is several orders 
of magnitude lower in the MS model (see text). The arrow shows the radius of a mass shell 
down to which the base of convective envelope penetrates at the end of FDU. 



